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f—N . We consider the minimal supersymmetric 5'O(10) model, where only one 10 

"^ \ and one 126 Higgs multiplets have Yukawa couplings with matter multiplets. This 

Q_i' model has the high predictive power for the Yukawa coupling matrices consistent 

ji. . with the experimental data of the charged fermion mass matrices, and all the Yukawa 

D i coupling matrices are completely determined once a few parameters in the model 

1^ ' are fixed. This feature is essential for definite predictions to the proton decay 

^ ■ rate through the dimension five operators. Althogh it is not completely general, 

k> , we analyze the proton decay rate for the dominant decay modes p -^ K^V by 

^ \ including as many free parameters as possible and varying them. There are two free 

parameters in the Yukawa sector, while three in the Higgsino sector. It is found that 

an allowed region exists when the free parameters in the Higgs sector are tuned so as 

to cancel the proton decay amplitude. The resultant proton lifetime is proportional 

to l/tan^/3 and the allowed region eventually disappears as tan/3 becomes large. 
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1 Introduction 

One particularly attractive idea for the physics beyond the standard model (SM) is the 
possible appearance of a grand unified theory (GUT) jTj, where the standard model gauge 
interactions are all unified into a simple gauge group. The successful gauge coupling 
unification of the minimal supersymmetric (SUSY) standard model (MSSM) at a scale 
Mg ^ 2 X 10^^ [GeV] j2] strongly suggests the circumstantial evidence of both ideas of 
SUSY and GUT, namely the SUSY GUTs. 

The most characteristic prediction of the (SUSY) GUTs is the proton decay. Normally 
in SUSY GUTs the proton decay process through the dimension five operators involving 
MSSM matters and the color triplet Higgsino turns out to be the dominant decay modes, 
since the process is suppressed by only a power of the Higgsino mass scale. Experimental 
lower bound on the proton decay modes p —>■ K^V through the dimension five operators 
is given by SuperKamiokande (SuperK) j3j, 

t{j) -^ K^V) > 2.2 X 10^^ [years]. (1) 

This is one of the most stringent constraints in construction of realistic SUSY GUT 
models. In fact, the minimal SUSY SU{5) model has been argued to be excluded from 
the experimental bound together with the requirement of the success of the three gauge 
coupling unification |lj [Sj. However note that the minimal SU{5) model predictions con- 
tradicts against the realistic charged fermion mass spectrum, and thus, strictly speaking, 
the model is ruled out from the beginning. Obviously some extensions of the fiavor struc- 
tures in the model is necessary to accommodate the realistic fermion mass spectrum. On 
the other hand, knowledge of the fiavor structure is essential in order to give definite 
predictions about the proton decay processes through the dimension five operators. Some 
models in which fiavor structures are extended have been found to be consistent with the 
experiments [7]. 

Recently, masses and fiavor mixings of neutrinos have been confirmed through the 
neutrino oscillation phenomena. This is the evidence of new physics beyond the standard 
model. One of the most promising candidates for new physics naturally accommodating 
the neutrino physics is the SUSY GUTs based on the gauge group 5*0(10). This is because 
in SO{10) models all the matters in the standard model together with additional right- 
handed neutrinos are unified into a single representation 16. Also the models can naturally 
explain the tiny neutrino masses compared to the electroweak scale through the seesaw 
mechanism p'.. Lots of models based on 5*0(10) have been proposed and extensively 
discussed. To be realistic, any models must accommodate all the experimental data of 
the fermion mass matrices and also be consistent with the experimental lower bound of 
the proton lifetime. As discussed above, knowledge of the fiavor structure in a given model 
is necessary in order to give a definite results on the proton lifetime. 

In this paper, we consider the so-called minimal SUSY 50(10) model proposed in P 
and analyzed in detail in fO] ^Ij ^2j, and we perform detailed analysis on the proton 



decay rate. The important fact is that this model has the high predictive power for 
the Yukawa matrices with leaving a few parameters free. Fixing these parameters, we 
can completely determine all the elements in the Yukawa matrices, whose information is 
essential for the definite predictions to the proton decay process through the dimension 
five operators. ^ Information of the Higgs sector in the model is also necessary for the 
definite predictions, and we include it in our analysis as free parameters. For simplicity, 
we assume all the soft masses are flavor diagonal. 

This paper is organized as follows: in the next section, we consider the minimal SUSY 
5*0 (10) model, and show the fact that the model has the high predictive power for the 
Yukawa matrices. In section 3 we analyze the proton decay processes with some free 
parameters and the predicted Yukawa matrices in the model. We show the distributions 
of the predicted proton lifetime for typical values of the free parameters. Section 4 is 
devoted to conclusions and discussions. 

2 Minimal SO (10) model and its predictions 

We begin with a review of the minimal SUSY 5*0(10) model proposed in P] and recently 
analyzed in detail in Ref. ^Oj ^I] ^2]- ^ Even when we concentrate our discussion on 
the issue how to reproduce the realistic fermion mass matrices in the 5*0(10) model, there 
are lots of possibilities for introduction of Higgs multiplets. The minimal super symmetric 
50(10) model is the one where only one 10 and one 126 Higgs multiplets have Yukawa 
couplings (superpotential) with 16 matter multiplets such as 



Wy = F/o'ieaOHie, + F/2'6l6a26^16, , (2) 



where 16j is the matter multiplet of the i-th generation, IOh and 126^;^ are the Higgs 
multiplet of 10 and 126 representations under 50(10), respectively. Note that, by virtue 
of the gauge symmetry, the Yukawa couplings, Yio and F126, are complex symmetric 
3x3 matrices. We assume some appropriate Higgs multiplets, whose vacuum expectation 
values (VEVs) correctly break the 50(10) GUT gauge symmetry into the standard model 
one at the GUT scale, Mg ~ 2 x W^[GeY]. Suppose the Pati-Salam subgroup ill], 
O422 = 5f/(4)c X SU{2)l X SU{2)ji, at the intermediate breaking stage. Under this 
symmetry, the above Higgs multiplets are decomposed as 10 -^ (6, 1, 1) + (1, 2, 2) and 
126 ^ (6,1,1) + (10,3,1) + (10,1,3) + (15,2,2), while 16 ^ (4,2,1) + (4,1,2). 
Breaking down to the standard model gauge group, 5f/(4)c x SU{2)ji -^ SU{3)c x f/(l)y, 
is accomplished by non-zero VEV of the (10, 1,3) Higgs multiplet. Note that Majorana 
masses for the right-handed neutrinos are also generated by this VEV through the Yukawa 



^The importance of flavor physics in dimension six proton decay processes is discussed in |13| . 

^In this section, we repeat almost the same discussions in Ref. ^U] (but input values we use are 
different) for the convenience of readers. If a reader is familiar with discussions here, the reader can skip 
to the the last paragraph of this section. 



coupling Yi26 in Eq. (0). In general, the SU{2)l triplet Higgs in (10,3, 1) C 126 would 
obtain the VEV induced through the electroweak symmetry breaking and may play a 
crucial role of the light Majorana neutrino mass matrix. This model is called the type II 
seesaw model, and we include this possibility in our model. 

After the symmetry breaking, we find two pair of Higgs doublets in the same repre- 
sentation as the pair in the MSSM. One pair comes from (1,2,2) C 10 and the other 
comes from (15, 2, 2) C 126. Using these two pairs of the Higgs doublets, the Yukawa 
couplings of Eq. (jSj) are rewritten as 

+ L, {Yli, vt) L, + (N^), (r^6 vr) (iV^). , (3) 

where If^, D'^, N'^ and E'^ are the right-handed SU(2)l singlet quark and lepton super- 
fields, Q and L are the left-handed SU{2)l doublet quark and lepton superfields, if "q ^'^d 
Hi2& are up-type and down-type Higgs doublet superfields originated from 10 and 126, 
respectively, and the last two terms is the Majorana mass term of the left-handed and the 
right-handed neutrinos, respectively, developed by the VEV of the (10, 3, 1) Higgs {vt) 
and the (10, 1,3) Higgs {vr). The factor —3 in the lepton sector is the Clebsch-Gordan 
coefficient. 

In order to keep the successful gauge coupling unification, suppose that one pair of 
Higgs doublets given by a linear combination H^q and -^"26 is light while the other pair 
is heavy (~ Mq)- The light Higgs doublets are identified as the MSSM Higgs doublets 
{Hu and Hd) and given by 

Hu = OIuHIq + PuHi2Q, 

Hd = adHt, + MU, (4) 

where au,d and (3u,d denote elements of the unitary matrix which rotate the fiavor basis 
in the original model into the (SUSY) mass eigenstates. Omitting the heavy Higgs mass 
eigenstates, the low energy superpotential is described by only the light Higgs doublets 
Hu and Hd such that 

Wy = {U% [a^Yll + m^) Hu Q, + {p% (a'Y'^^ + /J'^F/^) Hd Q, 

+ L, {Yli, vt) L, + {N% {Yli,vn) (iV^), , (5) 

where the formulas of the inverse unitary transformation of Eq. (0}, H^^ = a'^''^Hu,d+ ■ ■ ■ 
and Hi2e = l3^''^Hu,d + ■ • •, have been used. 



Providing the Higgs VEVs, Hu = f sin/3 and H^ = v cos (3 with v ~ 174.1[GeV], the 
quark and lepton mass matrices can be read off as 

Mu = CioMio + Ci26Mi26, 

Md = CioMio - 3Ci26Mi26, 
Me = Mio-3Mi26, 

Mt = ctMuq, 

Mr = CrMuG, (6) 

where M„, M^, M^, M^, Mt and Mr denote up-type quark, down-type quark, neutrino 
Dirac, charged-lepton, left-handed neutrino Majorana, and right-handed neutrino Majo- 
rana mass matrices, respectively. Note that all the quark and lepton mass matrices are 
characterized by only two basic mass matrices, Miq and Muq, and four complex coeffi- 
cients cio, Ci26, ct and cr, which are defined as Mio = Yiq a'^v cos j3, M126 = ^126 P'^'v cos /?, 
Cio = (a'"/a'^) tan/3, cue = (/?"//3'^) tan/?, ct = vt / {P'^v cos f3) , and cr = vr/{(3'^vcosP), 
respectively. 

The mass matrix formulas in Eq. (jH)) leads to the GUT relation among the quark and 
lepton mass matrices. 

Me = Cd {Md + kMu) , (7) 

where 

3Cio + Ci26 

Cd = , 

ClO — Ci26 

4 



3Cio + Ci26 

Without loss of generality, we can begin with the basis where Mu is real and diagonal, 
Mu = Du- Since Md is the symmetric matrix, it can be described as Md = V^^m ^d K^km 
by using the CKM matrix Vckm and the real diagonal mass matrix Dd- ^ Considering 
the basis-independent quantities, tr [MJMg], tr[(MJMe)^] and det[MJ'Me], and eliminating 
\cd\, we obtain two independent equations, 

tr[Me^Me] y ^ tr[(Me^Me)^] .g. 

ml + m?^ + ra^ j mj + m^ -|- m^ ' 

tr[Me^Me] V ^ det [Me^Me] . . 

9 I 9 I 9 I 9 9 9' V "J 

m^ + I^Ji + ^T / '^e "^^ "^T 



®In general, Md — ^ckm ^d ^ckm tiy using a general unitary matrix T^ckm = 
g2agi/3T3g27T8y^j^j^gi/3 73^17 Tg ^ -y^g omit the diagonal phases to keep the free parameters in the model 
minimum as possible. We can check that the final results of our analysis for the proton lifetime remains 
almost the same even if the diagonal phases are included. 



where Mg = VQ^^DdVQ^^ + kD^- With input data of six quark masses, three angles 
and one CP-phase in the CKM matrix and three charged lepton masses, we can solve 
the above equations and determine k and |q|, but one parameter, the phase of q, is left 
undetermined [UlI. The original basic mass matrices, Mio and M126, are described by 

MlO = ^^^FcKM^.T^^KM + ^^^., (11) 

Ml26 = ^^^FcKM^.^CKM-^^^.- (12) 

Here the unitary matrix I^ckm is reltated with the conventional Vckm as 

FcKM = e^'-e'^^^e'^'^WcKMe'^'^'e'^'^^ (13) 

However, in this paper we adopt these phases a, /?, 7, /?' and 7' to zero or tt. So Mio and 
Mi26 are the functions of a, the phase of q, with the solutions, \cd\ and k, determined by 
the GUT relation. 

Now let us solve the GUT relation and determine |q| and k. Since the GUT relation 
of Eq. ((7j) is valid only at the GUT scale, we first evolve the data at the weak scale 
to the ones at the GUT scale with given tan/5 according to the renormalization group 
equations (RGEs) and use them as input data at the GUT scale. Note that it is non- 
trivial to find the solution of the GUT relation, since the number of the free parameters 
(fourteen) is almost the same as the number of inputs (thirteen) . The solution of the GUT 
relation exists, only if we take appropriate input parameters. Therefore, in the following 
analysis, we vary two input parameters, ttIs and CP-phase 5 in the CKM matrix, within 
the experimental errors so as to find the solution. We take input fermion masses at Mz 
as follows (in GeV): 

ruu = 0.00233 , m^ = 0.677 , m^ = 176, 

rud = 0.00469 , m^ = 3.00, 

rUe = 0.000487 , m^ = 0.103 , m^ = 1.75. 

Here the experimental values extrapolated from low energies to Mz were used jjl^, and 
we choose the signs of the input fermion masses as — for rric, rrid and m^, and -|- for m„, 
irtt and r/if,. For the CKM mixing angles in the "standard" parameterization, we input 
the center values measured by experiments as follows [H 



_ 0.219 + 0.226 _ 0.037 + 0.043 _ 0.002 + 0.005 

■512 — , S23 — , Si3 — . 

In the following, we show our analysis in detail in two cases tan/3 = 2.5 and 10. The lower 
bound on tan/5 > 2.4 comes from the LEP Higgs searches J7]. We take m^ = 70.1 [MeV] 
at Mz scale for both cases of tan/5 = 2.5 and 10, but take 5 = 132° and 6 = 108° for 



tan/3 = 2.5 and 10, respectively. After RGE evolutions to the GUT scale, fermion masses 
(up to sign) and the CKM matrix are found as follows (in GeV): For tan/3 = 2.5, 

TTiu = 0.00126 , TTic = 0.365 , m* = 184, 
rrid = 0.00134 , m, = 0.0200 , nib = 1-07, 



and 

Vckm(Mg) 

and for tan/5 = 10 



and 



^ckm(Mg) 



me = 0.000323 , m^ = 0.0682 , m^ = 1.16, 



/ 0.975 0.222 -0.00188 - 0.00208i 

-0.222 - O.OOOlOli 0.974 + 0.000142? 0.0321 

V 0.00896 - 0.00203i -0.0308 - 0.000468i 0.999 



m„ = 0.000980 , m^ = 0.285 , m^ = 113, 
rrid = 0.00135 , m, = 0.0201 , nib = 0.996, 
me = 0.000326, m^ = 0.0687 , m^ = 1.17, 



/ 0.975 0.222 -0.000940 - 0.00289i 

-0.222 - 0.000129i 0.974 + 0.000124? 0.0347 

V 0.00864 - 0.00282i -0.0337 - 0.000647i 0.999 



in the standard parameterization. By using these outputs at the GUT scale as input 
parameters, we solve Eqs. (jH)) and (fTO|l . The contours of solutions of each equations are 
depicted in Figure^for tan/3 = 2.5. The crossing points of two contours are the solutions. 
As an example, we list a solution (at the upper crossing point) 



for tan/5 = 2.5, and 



-0.00675 + 0.000309Z , 

\cd\ = 5.99 , (14) 



-0.0103 + 0.000606? , 

\cd\ = 6.32 , (15) 



for tan /3 = 10. 

Once the parameters, \cd\ and k, are determined, we can describe all the fermion mass 
matrices as a functions of a from the mass matrix formulas of Eqs. ©, (fTT|) and (fT^ . 
Interestingly, in the minimal 5*0(10) model even light Majorana neutrino mass matrix, 
Mj,, can be determined as a function of a, ct and cr through the seesaw mechanism 



Ml, = —M^^M^^Md + Mt . The case where the the first term dominates is called type 
I seesaw model, while the case where the the second term dominates is called type II 
seesaw model. Each cases have been analyzed in detail in [10] and [TT], respectively, and 
discussed the consistency with the current neutrino oscillation data. Recently the general 
case has been analyzed, it is found that the minimal 5*0 (10) model is not so good to fit all 
the neutrino oscillation data J2|- This means that the structure of the model is somewhat 
too restrictive for the neutrino sector. It would be inevitable to extend the model at least 
for the neutrino sector in order to make the data fitting much better. Note that there 
are lots of possible ways to minimally extend the model only for the neutrino sector but 
keep the predictive power for the charged fermion mass matrices. Remember that there 
are two free parameters relevant to the neutrino sector, vt and vr-, in the model. If vt is 
smaller than the mass scale of the neutrino oscillation data, and if vr is the GUT scale^ 
the resultant light neutrino mass eigenvalues through the seesaw mechanism are too small 
to be compatible with the scale of the neutrino oscillation data. In such a case, we have 
to extend the neutrino sector. For instance, we can introduce an additional 126 multiplet 
which admits to obtain the VEV only in the SU{2)i triplet direction, w^. Now suppose 
that the type II seesaw mechanism works and we obtain the light neutrino Majorana 
mass matrix as Y(2Qv'rp, where Y/je ^^ newly introduced Yukawa coupling. If each elements 
in y/ge ^^^ much smaller than that in Yi2q, our analysis above remains correct and the 
predictive power for the charged fermion mass matrices is maintained. In the following 
analysis, we assume such a minimal extension of the model. 

3 Proton decay via dimension five operator 

The Yukawa interactions of the MSSM matter with the color triplet Higgs induces the 
following Baryon and Lepton number violating dimension five operator 

W = Cf^Q'Q^Q^LK (16) 

Here the coefficients are given by the products of the Yukawa coupling matrices and the 
(effective) color triplet Higgsino mass matrix, and are model dependent. In the minimal 
50(10) model, the coefficients are given by the products of two basic Yukawa coupling 
matrices, Y\q and Yviio^ and the effective 2x2 color triplet Higgsino mass matrix, M^, 
such as [IH] 

As discussed in the previous section the Yukawa coupling matrices, Fio and 1^26, are 
related to the corresponding mass matrices Mio and Mi26 such that 



'This is natural if the S0{1Q) group is broken down to the standard model one at the GUT scale. 



1^126 = ^/''' . Mi26, (18) 

jj^v sm p 

with f ~ 174.1 [GeV]. Here a" and /9" are the Higgs doublet mixing parameters introduced 
in the previous section, which are restricted in the range |a"p + |/9"P < 1. Although these 
parameters are irrelevant to fit the low energy experimental data of the fermion mass 
matrices, there are theoretical lower bound on them in order for the resultant Yukawa 
coupling constant not to exceed the perturbative regime. Since Cio, Cue, ^lo and M126 
are the functions of only a, we can completely determine the Yukawa coupling matrices 
once 0", a" and /3" are fixed. In order to obtain the most conservative values of the proton 
decay rate, we make a choice of the Yukawa coupling matrices as small as possible. In the 
following analysis, we restrict the region of the parameters in the range (a")^ + (/3")^ = 1 
(we assume a" and /S** real for simplicity). Here we present examples of the Yukawa 
coupling matrices with fixed a = n. For tan/? = 2.5 with a" = 0.031, we find 

0.000839 + 2.79 x 10"^ i 0.00151 - 0.0000265 i 0.000692 - 0.000818 i 
Yio=\ 0.00151- 0.0000265 i 0.00479 + 0.0000811 i -0.0128 + 3.17 x 10"^^ 
0.000692 - 0.000818 i -0.0128 + 3.17 x 10"^ i 0.525 - 0.0420 i 

(19) 
-0.0000613 - 2.17 X 10"^ i -0.000111 + 1.94 x 10"^ i -0.0000508 + 0.0000600 i 
Fi^R = I -0.000111 + 1.94 X lO^S -0.000428 - 2.46 X lO^S 0.000941 - 2.33 x 10"^ i 
-0.0000508 + 0.0000600 i 0.000941 - 2.33 x 10"^ i 1.42 x 10"^ + 0.00132 i 

(20) 
and for tan/? = 10 with a" = 0.111, 

0.00101 + 1.87 X 10-6 i 0.00179- 0.0000439 i 0.000348 - 0.00125 i 
rio=| 0.00179- 0.0000439 i 0.00541 + 0.00014H -0.0154 + 5.21 x 10"^^ 
0.000348 - 0.00125 i -0.0154 + 5.21 x 10"^ i 0.530 - 0.0567i 

(21) 
-0.000244-5.21 X 10-^ i -0.000436 + 0.0000107? -0.0000847 + 0.000305 i 
Yi26=\ -0.000436 + 0.0000107i -0.00164 - 0.0000154i 0.00375 - 1.27 x 10"^ 
-0.0000847 + 0.000305 i 0.00375 - 1.27 x 10"^ i -0.000684 + 0.00626 i 

(22) 
For the effective color triplet Higgsino mass matrix, we assume the eigenvalues being 
the GUT scale, Mq = 2 x 10^^ [GeV], which is necessary to keep the successful gauge 
coupling unification. Then, in general, we can parameterize the 2x2 mass matrix as 

Mc = Mch X f/, (23) 

with the unitary matrix, 

JJ^^^^aJ COS0 smfii \ ^^Vg _ (24) 

\ — sm U cos U ^ 



Here we omit an over all phase since it is irrelevant to calculations of the proton decay 
rate. Now there are five free parameters in total involved in the coefficient C^l , namely, 
0", a„, 6*, if and ip' . Once these parameters are fixed, C^ is completely determined. 

The proton decay mode via the dimension five operator in Eq. (fTT)|l with the Wino 
dressing diagram is found to be dominant, and leads to the proton decay process, p -^ 
K^V. The decay rate for this process is approximately estimated as (in the leading order 
of the Cabibbo angle A ~ 0.22) 

X 5.0 X 10^^ [years" VGeV]. (25) 

Here the first term denotes the phase factor and the hadronic factor, (3h = 0.0096 [GeV^'^] 
given by lattice calculations \^. A^ ~ 0.32, As — 0.93 are the long-distance and the 
short-distance renormalization factors about the coefficient C^ , respectively. ^° The 
third term in the first line comes from the Wino dressing diagram, and ms is a typical 
sparticle mass scale multiplied by the ratio of a sfermion and Wino. In the case with the 
mass hierarchy between the sfermions and the Wino (m/ ^ M2), we find ms ~ mj x 
{mf/M2). In the following numerical analysis, we take rfif = 1 [TeV] and M2 = 100 [GeV]. 
Now we perform numerical analysis. Note that because of the very constrained flavor 
structure of the minimal SUSY 5*0(10) model we can give definite predictions for the 
proton decay rate once the five parameters in the above are fixed. For a specific choice 
of the Yukawa coupling matrices in the minimal S'0( 10) model with the type II seesaw, 
the proton decay rate has been calculated in |2I]. In our analysis, we make no such a 
specific choice, and perform detailed analysis in general situations of the minimal 5*0(10) 
model by varying the above five free parameters. The result for tan/3 = 2.5 is presented 
in Figure |2l Here the distributions of the proton lifetime (log years) for arbitrary choices 
of the five free parameters (normalized by 1) is depicted. We can see that some special 
sets of the free parameters can result the proton lifetime consistent with SuperK results. 
In that region, cancellation in the second line in Eq. ()25|) occurs by tuning of the free 
parameters in the Higgsino mass matrix. Note that number of free parameters is not 
enough to cancel both of the process p -^ K~^T^ (dominant mode) and p — > K^V^ (sub- 
dominant mode), and thus the proton lifetime has an upper bound in the 5*0(10) model. 
For tan/? = 10, we obtain the same figure depicted in Figure El but the lifetime is scaled 
by roughly (2.5/10)^, which is consistent with the naive expectation that the lifetime is 
proportional to 1/ ia-v? (3. Whole region is excluded in the case with tan/5 = 10. 

^°As suggested in Ref. [201 1 it might be proper to use the renormahzation factors A^, As in [201, which 
directly treats the renormahzation of the Wilson coefficients itself. But here, we adopt the use of the 
conventional factors Al, As to compare our results to the previous ones. 
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In the case of nondegenerate masses of Mc, the parameters increase from three to 
5 + 1 (the last 1 is the ratio of masses). However, the results only slide by the square of 
this mass ratio and do not show the special cancellation. 

4 Conclusion and discussions 

We have discussed the minimal SUSY SO (10) model, which can reproduce the realistic 
charged fermion mass matrices with only one parameter left free. This model has high 
predictive power for the fermion Yukawa coupling matrices, and they are completely 
determined once a few parameters in the model fixed. This feature is essential for definite 
predictions of the proton decay rate vis the dimension five operators. Including additional 
3 free parameters in the effective Higgsino mass matrix with mass eigenvalue being the 
GUT scale, we have analyzed the proton decay rate by varying the five free parameters 
in total. We have found that for tan/5 = 2.5 some special sets of the parameters predicts 
the proton decay rate consistent with the SuperK results, where the cancellation for 
the dominant modes of the proton decay amplitude occurs by tuning of the parameters. 
Although there exists the allowed region, it is very narrow. Our result is consistent with 
the one in the previous work [2Tj for only one specific choice of the Yukawa coupling 
matrices. It has been found that the resultant proton decay rate is proportional to tan^ (3 
as expected and the allowed region eventually disappears as tan/? becomes large, even for 
tan/3 = 10. 

There are some theoretically possible ways to extend the proton lifetime. One way 
is to adopt a large mass hierarchy between the sfermions and the Wino as can be seen 
in Eq. (j^ . The proton lifetime is pushed up according to the squared powers of the 
mass hierarchy, and the allowed region becomes wide. How large the hierarchy can be 
depends on the mechanism of the SUSY breaking and its mediation. When we assume the 
minimal supergravity scenario, the cosmologically allowed region j22] consistent with the 
recent WMAP satellite data [21] suggests that masses between sfermion masses the Wino 
is not so hierarchical and the value we have taken in our analysis seems to be reasonable. 
Another way is to abandon the assumption of Higgsino degeneracy at the GUT scale, and 
to make the mass eigenvalues of the effective colored Higgsino mass matrix heavy. We 
can examine this possibility based on a concrete Higgs sector. However this seems to be 
a very difficult task even if we introduce a minimal Higgs sector in the minimal 5*0(10) 
model discussed in jTHj [21], since there are lots of free parameters in the Higgs sector. 
Furthermore, even in the minimal Higgs sector, there are lots of Higgs multiplets involved 
and the beta function coefficients of the gauge couplings are huge. It seems to be very 
hard to succeed the gauge coupling unification before blowing up of the gauge couplings. 
Therefore, the assumption that all the Higgs multiplets are degenerate at the GUT scale 
would be natural. Consequently our results show the typical properties of SO(IO) GUT 
but are not exhaustive. Also there is possibility to vary GUT phases a, /3, 7, (3' and 7' in 
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Eq. (unD. 
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Figure 1: Contour plot on complex K-plane in case of tan/3 = 2.5. The vertical line and 
the circle correspond to the solutions of Eqs. ^ and (fTUI) . respectively. 
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Figure 2: The distributions of the proton lifetime (log years) for tan/3 = 2.5 in arbitrary 
five parameter choices (normalized by 1). The green region (right hand side) is consistent 
with experiment. 
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Figure 3: Same as in Figure |2l but for tan/? = 10. 
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